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Telephone Building Energy Consumption and Control (TELBECC) pro- 
gram has been developed to accurately and efficiently analyze environmental 
control and energy use in telephone company buildings. The program simulates 
various operational plans to determine the relative energy and cost savings. 
By analyzing the operation of the heating, ventilation, and air conditioning 
system as it regulates a changing environment, TELBECC calculates the 
heating and cooling load, dry-bulb temperature, and relative humidity in the 
building. The user specifies the building's dry-bulb temperature limits, which 
are the control variables for the program analysis. The simplified computa- 
tional procedure of the program incorporates a recursive scheme using time 
series to perform the necessary calculations. The results of the computations 
can be obtained for different periods: the quarter hour, hour, day, or month. 
Energy consumption and control in several equipment buildings located in 
three different geographical areas have been analyzed by TELBECC. Analysis 
and comparison of the resulting data demonstrate the advantages of the 
program. 

I. INTRODUCTION 

An ambitious energy-cost savings program has been instituted to 
reduce energy use in telephone company buildings. In recent years, 
telephone companies have saved energy mainly by redesigning and 
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retrofitting buildings to operate and maintain environmental control 
equipment at peak performance, to turn out unneeded lights, to reduce 
heating and cooling losses, and the like. Further energy and cost 
savings, although requiring additional capital investment, could be 
achieved through modification of environmental control systems, pur- 
chase of sophisticated microprocessors for more efficient control of 
building Heating, Ventilating, and Air conditioning (HVAC) systems, 
and installation of alternative energy sources, such as solar power and 
wind power. Before adopting any conservation plans that require 
appreciable capital investment, however, we should make a thorough 
economic evaluation. Such an evaluation can be carried out by corre- 
lating the changing operating characteristics of a building with selected 
energy conservation plans. This procedure would enable us to pinpoint 
the most economical operating strategies. 

There are several commercially available computer programs to 
perform this type of analysis, such as DOE II, ESP, and BLAST; 1 - 2 
most, however, are proprietary. These complex programs can analyze 
any of a broad spectrum of commercial, industrial, and residential 
buildings. But, because of their versatility, they require large computer 
systems, extensive data preparation, and high costs. The use of energy 
in the majority of telephone company equipment buildings, which are 
small, single-story structures varying in area from 1500 to 10,000 
square feet, can be best evaluated by a more focussed computer 
program. 

This paper describes a new computational method and computer 
program called Telephone Building Energy Consumption and Control 
(TELBECC). This program simulates building operations and quickly 
evaluates numerous energy conservation plans and cost-saving strat- 
egies under variable weather conditions [according to standard hourly 
Test Reference Year (TRY) weather data 3 ]. The program can evaluate 
energy consumption for intermittent or proportional HVAC plant 
operation, economy cycle operation, enthalpy cycle operation, and 
wideband temperature operation with no heating or cooling between 
preset room temperature limits. Also, the program can calculate the 
optimum building orientation and U factor (heat transmission char- 
acteristics) of the outside walls and roof, chiller and heater plant size, 
dry-bulb temperature and Relative Humidity (RH) variations, and 
quantity of water required to maintain 20-percent minimum RH 
during economy cycle operations. 

II. PROGRAM DESCRIPTION 

We can derive the heating or cooling load in an enclosed building 
from the following considerations: 
1. Conduction of heat through the building walls and roof. 
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2. Permeation of outside air through the building envelope. 

3. Internal heat generation from equipment, lights, and people. 

4. Direct solar radiation through windows and skylights (fenestra- 
tion). 

Since most operating company equipment buildings have few win- 
dows, the program does not consider item 4 in the present version. 

A Constant-Air- Volume (CAV) supply fan system typically controls 
the air temperature of a building. The building engineer normally sizes 
the fan system using the elementary steady-state heat balance, which 
takes into account the internal heat loads, outside air temperatures, 
and solar radiation in conjunction with the U factor of the building 
envelope. In general, this conservative approach produces fan systems 
that are oversized and therefore inefficient. To find a smaller, and 
perhaps more efficient, fan capacity, a TELBECC user selects different 
fan capacities for analysis by the program. The program generates 
data on the space temperatures, relative humidity, peak heating and 
chiller loads, and the hours of system operation for the different fan 
capacities that can be used to find an optimum air supply fan system. 

For comfort, a limit is imposed on the difference between HVAC 
supply and return air temperatures. For cooling, this temperature 
difference is —20 degrees F; for heating, +40 degrees F. These default 
values may be overridden by the user. With an environmental dry- 
bulb temperature standard specified, the program computes the re- 
quired operation of the HVAC. 

The user can specify one of two basic ways to operate the HVAC: 
intermittent operation or proportional control. With intermittent 
operation, the HVAC does not supply any heating or cooling when the 
dry-bulb air temperature is within the wideband temperature range. 
Reaching or exceeding either wideband temperature limit activates 
the HVAC. The HVAC stays on and does not deactivate until the dry- 
bulb air temperature reaches 3 degrees F above the lower limit of the 
wideband temperature range for heating and 3 degrees F below the 
upper limit of the wideband temperature range for cooling. The TEL- 
BECC user can reset the numerical values of the throttling range if a 
different range is appropriate. The proportional control plan operates 
by continuously adjusting the supply and return air temperature 
difference in increments of 1 degree F to satisfy the instantaneous 
building heating or cooling load. This plan follows the building load 
to closely track the lower and upper limits of the wideband temperature 
range with essentially no throttling. When selecting a dual or extended 
wideband temperature standard (that is, one with different wideband 
limits for occupied and unoccupied times), the HVAC activates before 
occupancy in order to reach the preset temperature standard. 

TELBECC calculates the heat added or removed by the HVAC 
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system in controlling the dry-bulb air temperature every quarter hour. 
In particular, when cooling is required, the sensible and latent loads 
on the chiller plant are simultaneously computed by incorporating any 
of three standard methods of fan system operation: 

1. Conventional operation, which is chiller operation with no econ- 
omizer. 

2. Chiller operation with a dry-bulb economy cycle. 

3. Enthalpy cycle. 

In conventional operation, the minimum quantity of outside air needed 
for ventilation is circulated. This mode is also used as a benchmark 
for the program. The dry-bulb economy cycle uses outside air for 
cooling whenever the outside dry-bulb air temperature falls below the 
maximum value. The default value is 55 degrees F, but the user can 
reset the value. The enthalpy cycle checks the enthalpy of the inside 
air and the outside air. If the outside air enthalpy is lower, 100 percent 
outside air is circulated to reduce the load on the chiller regardless of 
the relative humidity. Otherwise, only the minimum quantity of out- 
side air required for ventilation is circulated. 

System control is based on dry-bulb air temperature and is not 
predicated on maintaining a particular value of relative humidity. 
Nevertheless, the program computes changes in relative humidity for 
the three methods of fan-system operation discussed above. The 
program summarizes the variation in relative humidity for the time 
period chosen by giving the number of hours the relative humidity is 
less than 10 percent, between 10 and 15 percent, between 15 and 20 
percent, between 20 and 55 percent, between 55 and 60 percent, and 
greater than 60 percent. 

In addition, since dry-bulb economy cycle operation generally calls 
for bringing in winter air with low humidity, the program calculates 
the quantity of water required for humidification. The operating 
company minimum standard of 20 percent RH in the inside air for 
dry-bulb economy cycle operation in winter is the basis for calculating 
the amount of water added to the air. 

III. TRANSIENT HEAT CONDUCTION THROUGH THE BUILDING 
ENVELOPE 

Weather conditions influence the heating and cooling load of a 
building by heat conduction through the structural and decorative 
materials of the exterior walls and roof, as well as by permeation of 
outside air and direct absorption of solar radiation through window 
areas. Since, as previously mentioned, most operating company equip- 
ment buildings have few windows, only heat conduction and permea- 
tion are treated in the computer program. The program must account 
for the heat storage effects of the structure, as well as the daily and 
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seasonal variation of the outside air temperature and solar radiation. 
We can account for these influences on the building by considering 
the building envelope elements as one-dimensional flat slabs or plates. 
We then obtain a solution to a partial differential equation with time- 
dependent boundary conditions. A classical analytical solution of this 
equation 4 produces a set of equations that require an inordinate 
quantity of computational effort and time, rendering the whole idea 
of performing the analysis impractical and uneconomical. However, 
the analytical solution can be recast into a simpler, though effective, 
computational scheme with a method first introduced by Mitalas and 
Stephenson, 5 which is ideally suited to calculation by computer. 

The inside-wall and roof-surface temperature TbeU) and the air 
temperature of the building T a (t), which are dependent on time, 
determine the environmental load due to convection.* T B e(0 is rep- 
resented in the form of the following time series: 



m-1 



T BE( = - 2 fcTBE,.,. + 2 OiTo^ + 2 aiT Bl _ i _ lt (1) 

i=i i=i ;=i 

where 

t = current time, 
Tbe, = inside-wall temperature of the building at the current time, 

7be,_, = inside-wall temperature of the building i time units prior 

to t, 
To,., = outside sol-air temperature 6 i time units prior to t, and 

T a/ _. = air temperture of the building i time units prior to t. 

For hourly temperature calculations, the number of terms, m, will 
rarely exceed 5, and for quarter-hour calculations, m will generally be 
less than 15. The recursive properties of the calculation make it 
extremely efficient and economical, especially when the operating 
characteristics of the building may need to be tracked for an entire 
calendar year. The coefficients 6„ a„ and a/ in eq. (1) are determined 
from the thermophysical properties of the structure. Only six values 
are needed to uniquely specify these coefficients: wall thickness, wall 
U factor, wall-weight density, effective heat-transfer coefficient of the 
inside- and outside-wall surface, and the time interval between 
successive calculations. The appendix presents the mathematical pro- 
cedure to evaluate these coefficients. Figures 1 and 2 show the math- 
ematical and physical models for deriving the coefficients. 
We can validate this simplified computational approach by compar- 



* Radiative interchange between inside building-wall surfaces is not included in the 
present version of the program. 
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Fig. 1 — Discrete and continuous functions. 
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Fig. 2 — Homogeneous flat slab. 

ing it with an exact solution given in the literature. 7 We can see this 
in Figs. 3 and 4 for several values of the inside-wall dimensionless 
convective heat-transfer parameter B\ ■ h\L/k and two limiting values 
of the outside-wall convective heat-transfer parameter Bo ■ hoL/k. In 
Fig. 3 the outside-wall convective heat transfer parameter Bo = 0; i.e., 
the surface x = L is insulated. In Fig. 4 the solution corresponds to B 
approaching 00; i.e., the surface x = L is maintained at a constant 
temperature. The initial and boundary conditions are indicated in the 
figures. The solid curves represent the exact condition given in Ref. 7, 
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Fig. 3 — Temperature response of front face of plate, < x < L, with back face < x 
< L, with back face x = L maintained at T = degree F after sudden exposure to 
uniform-temperature convective environment T a = 1 degree F at x = 0. Sampling 
interval A is one min. 



and the dots represent the numerical values computed by the time 
series in eq. (1). The sampling interval for this example is A = 1 min. 
or, in terms of dimensionless time, kA/L 2 = 0.001. We can see that 
after some time has elapsed the exact solution and the time -series 
solution match identically. For the problem considered here at t = 0, 
the ambient convected temperature T a is suddenly raised from T a = 
to T a = 1; i.e., the boundary condition is a step function. However, 
since the development of the time series assumes a linear variation 
between time intervals, as stated in the appendix, the solution resem- 
bles an initial ramp followed by a constant value, as shown in Fig. 5. 
Once the effect of the initial ramp input diminishes after about 10 
sampling intervals, the solution coincides with the exact solution. This 
characteristic of the time series is not a problem here, since instanta- 
neous changes of air temperatures inside and outside the building do 
not occur. 



IV. CALCULATION OF BUILDING AIR TEMPERATURE AND ENERGY 
USE 

The air temperature of the building is obtained through the following 
equation for the heat balance within the building: 

Qairit) = QequiptU) + 9lighta(0 + <7people(0 

+ <7innitration(t) + 9walls(0 + <?HVAcU), (2) 
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Fig. 4— Temperature response of front face of plate < x < L, with insulated back 
face x = L after sudden exposure to uniform-temperature convective environment T« = 
1 degree F at x = 0. Sampling interval A in one min. 
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Fig. 5— Convection temperature environment at inside boundary surface x = 0. 

where q aiT (t) (Btu/hr) represents the sensible thermal-energy convec- 
tion rate of the inside air, resulting in a change in the overall dry-bulb 
air temperature. The other terms on the right represent the rate at 
which heat is convected to the air from the following: equipment heat 
dissipation; lighting; people; inadvertent infiltration of outside air; the 
exterior walls, floor, and roof of the building; and the building's HVAC 
control system. 

To evaluate the air temperature T a (t) from eq. (2), each term is 
expressed as the difference in temperatures between the air and the 
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heat-convecting medium. We can determine T a (t) from the following 
differential equation: 

f™ ^JT- = H ELA E L[T EL (t) ~ TM] + pcm[T (t) - T e (t)] 
at 

+ H BE A BE [T BE (t) - T a (t)] + P cQ[DTsAt)), (3) 

where 

T a (t) = average building dry-bulb air temperature; 
TelU) = combined average temperature of the equipment, lights, 

and people; 
To(t) = outside dry-bulb air temperature; 
T BE {t) = inside surface temperature of the building envelope; 
DTsp(t) = difference between the air-supply temperature and the air 
temperature of the building; 
p = air density; 
c = air specific heat; 
v = volume of the building space; 
#el = heat-transfer coefficient between the air and equipment; 
-Ael = average surface area of equipment; 

m = rate at which outside air infiltrates the building; 
H BE = heat-transfer coefficient between the building envelope 

and air; 
A BE = surface area of the building envelope; and 
Q = air-supply rate of the HVAC system fan units. 

Only one of the temperatures in eq. (3) is presumed known: the out- 
side-air dry -bulb temperature, To(t). The other four tempera- 
tures— T a (0, T EL (t), T BK (t), and DT SP (t)— are coupled; therefore, ad- 
ditional equations are needed for their resolution. 

We can consolidate our terms into two other equations. When we 
combined the heat gain generated by the equipment, lighting, and 
people into a single term for heat dissipation per unit of building floor 
area W(t)(W /ft 2 ), one additional equation can be written as 



?el ^^ = 3.41 A t W(t) + H EL A EL [T a (t) - T Eh (t)], (4) 



where 



Cel = heat capacity of the equipment, lights, and people; 
A t = building floor area; and 
W(t) = combined heat dissipation of the equipment, lights, and peo- 
ple per unit of building floor area. 

A third equation coupling the building envelope temperature T BE (t) 
and the building air temperature T a (t) is needed. A likely equation 
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would be the time series given by eq. (1). But before this can be 
applied, both eqs. (3) and (4) must also be recast in the form of time 
series. This is easily done by using the properties of the z transform 
and the procedure already delineated in the apperidix. The time-series 
solution of eqs. (3) and (4) assumes the form 

T a , = s7V, + [ T \y j fo L T EL + toTo + r BE T BE + r SP DT SP ] t 

+ f 1 ~ S( ^2 + J [ t elT el + toTo + tbeTbe + r SP DT SP ] t - U (5) 
and 

T EL , = sTel,-, + (AT + Jt~ ^ [^ + *- T J* 

[l--s(Afj-J^] . ,. 

+ -75 It w VV + Taiajt-l, (o) 

where 

A = sampling time interval, 
tel = HelAel/pcv, 
to = m/u, 
tbe = HbeAbe/pcv, 

T SP = Q/V, 

T = TEL + TO + TBE + TSP, 

s = EXP(-At), 

T w = A{/CeU 

f a = HeL-^El/CeL, 

f = f w + t a , and 
s = EXP(-Af), 

and the subscripts t and t—1 indicate that the temperature is evaluated 
at times t = nA and t= (n— 1)A. 

The calculation of the humidity ratio in the space is also formulated 
in terms of a first-order linear differential equation in time similar to 
eqs. (3) and (4). This equation is also recast in the form of a time 
series [see eqs. (5) and (6)]. Knowing the humidity ratio and the dry- 
bulb air temperature, we can find the relative humidity by employing 
standard psychometric formulae. 6 

As previously noted in Section II, the system of eqs. (1), (5), and (6) 
permits controlling the dry-bulb air temperature in two basically 
different ways, intermittent operation and proportional control, and 
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determines the hours of operation. We use these equations for inter- 
mittent HVAC operation where we assume that the HVAC system 
holds the supply and return air-temperature difference constant at 
DT S p = —20 degrees F whenever cooling is required and at DT S p = 40 
degrees F whenever heating is required. The hours of operation, and 
hence the total quantity of heat removed or added to satisfy the 
imposed dry-bulb temperature standard, can then be determined. For 
the proportional control plan, eqs. (1), (5), and (6) are also used to 
calculate not only the hours of HVAC operation but also the numerical 
value of the supply and return air-temperature difference DTsp(t), 
which in general varies continuously for this mode of control. The 
variation in time of the numerical value of DT SP (t) is determined by 
just satisfying the instantaneous building environmental load. When 
the HVAC system is activated, the proportional control plan closely 
follows the lower limit (for heating) or the upper limit (for cooling) of 
the building wideband temperature range. 

Once DTsp(t) and the hours of operation are known, the program 
calculates the heat added or removed from the building by the HVAC 
system during every quarter hour and for whatever other period is of 
interest, e.g., monthly. As a corollary, we can estimate the environ- 
mental control system energy use, assuming the following HVAC 
system characteristics: 

1. For chiller operation, a constant Coefficient Of Performance 
(COP) supplied by the user, together with the quantity of heat removed 
from the building air, characterizes its energy requirements. 





Table 1- 


—Equipment buildings analyzed 






Energy Consumption and Control 




Wideband Temperature 

(°F) 


Control 




Case 


Occupied U 
Times 


noccupied 
Times 


Geographic Location 


1 

2 
3 
4 
5 


65-80 
65-80 
65-80 
65-80 
65-80 


65-80 
65-80 
60-85 
65-80 
65-80 


Intermittent 
Proportional 
Intermittent 
Intermittent 
Intermittent 


New York City 
New York City 
New York City 
New Orleans 
Phoenix 


Building Parameters 




Factor 






Parameter 



Size (L X W x H) 60 ft x 40 ft x 13 ft 

Average heat transmission U = 0.25 Btu/hr — ft 2 

Occupancy time 8 a.m. to 6 p.m. 

Fan support rate 7400 CFM 

Ventilation capacity 150 CFM 

Static fan pressure 2 in. of water 

Internal heat load 15 W/ft 2 

Economy cycle temperature limit <65 degrees F 
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2. For fan operation, fan power is calculated from the following 
equation: 



Fan power (kW) = 2.487 X 10 -4 X Q X AP, (7) 



where 



Q(CFM) = air flow delivered by the fan; and 

AP = static pressure head of the fan in inches of water. 

By multiplying the fan power by the total hours of fan operation, we 
can obtain the total energy use (kWh). 

3. Humidification costs are based on supplying energy at the rate 
of 1000 Btu per pound of water added to the supply air stream. Costs 
are derived from the unit cost of energy, such as electricity ($/kWh), 
natural gas ($/1000 ft 3 ), and fuel oil ($/gal), which is supplied by the 
user. An 80 percent efficiency rate is assumed for these energy sources. 

4. Heating costs are similarly calculated by the unit cost. An 80 
percent efficiency rate is also used in these calculations. 

V. ILLUSTRATIVE EXAMPLES OF ENERGY CONSUMPTION AND 
CONTROL 

Different geographic loctions of equipment buildings, dual or ex- 
tended wideband temperature limits, and the method of HVAC control 
(intermittent or porportional) are considered in several variations. 
Table I gives this information along with some of the more salient 
building parameters. The results of the calculation are summarized by 
month in Tables II through VI for the cases specified in Table I. We 
assume here that a conventional cooling system, consisting of a chiller 

Table Mb — Intermittent control of building space air temperature — 
New York City, case 1 in Table I 

Number of Hours at Specified Relative Humidity (No Humidity Control) 

<10% 10-15% 15-20% 20-55% 55-60% >€0% 

Conv (no econ) 697.00 1090.00 1086.75 5886.25 61) 0.0 

Economy 815.25 1048.75 1029.75 5866.25 0.0 0.0 

Enthalpy 697.50 1101.00 1144.00 5817.50 0.0 0.0 

Estimated Operating Cost for Cooling at $0.10/kWh for Electricity 
(Chiller COP = 3.50) 

Conv (no econ) $9687 for 86870 kWh (Fans = 18669 kWh, chiller = 68201 kWh) 
Economy $8622 for 56221 kWh (Fans = 18669 kWh, chiller = 37552 kWh) 

Enthalpy $7534 for 75342 kWh (Fans = 18669 kWh, chiller = 56673 kWh) 

Estimated Operating Cost for Humidification (20% min) and Heating at $0.10/kW"h for 

Electricity 

Humidification $1959 for 19580 kWh 
Heating $0 for kWh 

Notes: Min space temp occurred on day 300, max space temp occurred on day 213, 
max cooling load occurred on day 224. 
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and air-handling unit, provides cooling. By comparing the different 
cases, we find some interesting results. 

Case 1 differs from case 2 in that the HVAC is intermittently 
controlled in case 1, but proportionally controlled in case 2. We see 
from Tables lib and Illb, for example, that the maximum cooling 
loads for both cases occur close in time [August 12 (day 224) and July 
31 (day 212)]. However, the maximum cooling load of 13 tons for the 
proportional control plan (Table Ilia), which compares favorably with 
the load of 14.2 tons for the intermittent control plan (Table Ha), 
reduces the required size of the chiller plant by 9 percent. We would 
expect such a reduction from using a control sequence that follows the 
load closely and minimally overshoots the dry-bulb air temperature. 
Also, a control plan that matches the fan capacity to the load would 
compare favorably in energy use with on-off fan operation. 

We can see in Table Ilia (in the column labeled "Space Temp") 
that the proportional control plan regulates the temperature to within 
one-tenth of a degree of the wideband temperature limit for the entire 
calendar year. For the economy cycle operation, the yearly electrical 
use of case 1 in Table lib is 56,221 kWh, and that of case 2 in Table 
Illb is 63,614 kWh. The chiller energy consumption for case 1, 37,522 
kWh, is larger than that for use 2, 31,591 kWh. The proportional 
control plan, which modulates the air-supply temperature, requires 
the fan to run continuously at maximum power for the entire year. 
This maximum use of the fan creates larger overall energy require- 
ments in spite of lower chiller energy use. However, a variable-air- 
volume system that modulates the fan supply rate to match loads 
should decrease the required fan power and significantly reduce total 
energy use. 

Table Illb — Proportional control of building space air temperature — 
New York City, case 2 in Table I 

Number of Hours at Specified Relative Humidity (No Humidity Control) 

<10% 10-15% 15-20% 20-55% 55-60% >60% 

Conv (No Econ) 869.75 1090.75 1118.75 4393.25 1011.25 276.50 
Economy 911.25 1092.50 1110.50 4354.00 1015.00 276.75 

Enthalpy 869.75 1095.50 1134.50 4337.75 938.75 383.75 

Estimated Operating Cost for Cooling at $0.10/kWh for Electricity 
(Chiller COP = 3.50) 

Conv (no econ) $9741 for 97400 kWh (Fans = 32023 kWh, chiller = 65383 kWh) 
Economy $6361 for 63614 kWh (Fans = 32023 kWh, chiller = 31591 kWh 

Enthalpy $8930 for 88296 kWh (Fans = 32023 kWh, chiller = 56273 kWh) 

Estimated Operating Cost for Humidification (20% min) and Heating at $0.10/kWh for 

Electricity 

Humidification $6022 for 60219 kWh 
Heating $0 for kWh 

Notes: Min space temp occurred on day 100, max space temp occurred on day 10, 
max cooling load occurred on day 212. 
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Tables IVa and IVb display monthly energy use for case 3 of Table 
I. Differing from case 1, this plan imposes dual wideband temperature 
limits. The wideband temperature limits for unoccupied times increase 
to 60 degrees F and 85 degrees F. Tables IVa and IVb show that this 
simple change with economy cycle operation reduces annual cooling 
energy use by 7 percent, from 56,222 kWh to 52,379 kWh. We can 
attribute this saving mainly to the lower chiller energy requirement, 
from 37,552 kWh to 34,237 kWh, and to a lesser extent to the smaller 
fan energy requirements, from 18,669 kWh to 18,141 kWh. 

Tables V and VI show the results of the simulation for the buildings 
located in New Orleans and Phoenix, cases 4 and 5 of Table I. We can 
see that the maximum chiller load, 14.4 tons, is the same for these 
diverse locations. The total energy required in these buildings for 
economy cycle cooling is also nearly equal, 85,027 kWh and 85,123 
kWh. Since the cooling load includes both sensible and latent energy, 
we can surmise from the degree-day totals that the dominant load on 
the system for Phoenix is sensible heat, and on that for New Orleans 
latent heat. The distribution of the relative humidity and the costs for 
humidification shown in Tables Vb and VIb tends to support these 
observations. 

We can see the advantages in running the TELBECC program to 
compare different control plans. For example, although the intermit- 
tent control plan appears to use less energy overall, the proportional 
control plan actually reduces the size of the chiller plant by 9 percent 
in the same locale under similar conditions. The higher overall costs 
can be attributed to continuous operation of the fan, which, if operated 
to match the load, would consume much less power and make the 

Table IVb — Intermittent control of building space air temperature — 
New York City, case 3 in Table I 

Number of Hours at Specified Relative Humidity (No Humidity Control) 

<10% 10-15% 15-20% 20-55% 55-60% >60% 

Conv (No Econ) 857.50 1205.75 1213.00 5483.75 O0 0.0 

Economy 976.75 1115.75 1138.25 5529.25 0.0 0.0 

Enthalpy 857.50 1241.00 1212.00 5449.50 0.0 0.0 

Estimated Operating Cost for Cooling at $0.10/kWh for Electricity 
(Chiller COP = 3.50) 

Conv (no econ) $8425 for 84247 kWh (Fans = 18141 kWh, chiller = 66106 kWh) 
Economy $5238 for 52378 kWh (Fans = 18141 kWh, chiller = 34237 kWh) 

Enthalpy $7278 for 72776 kWh (Fans = 18141 kWh, chiller = 54635 kWh) 

Estimated Operating Cost for Humidification (20% min) and Heating at $0.10/kWh for 

Electricity 

Humidification $2380 for 23802 kWh 
Heating $0 for kWh 

Notes: Min space temp occurred on day 22, max space temp occurred on day 224, 
max cooling load occurred on day 224. 
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Table Vb— Intermittent control of building space air temperature- 
New Orleans, case 4 in Table I 

Number of Hours at Specified Relative Humidity (No Humidity Control) 

<10% 10-15% 15-20% 20-55% 55-60% >60% 

Conv (No Econ) 14.25 395.50 441.25 7909.00 o!b~ 0.0 

Economy 60.0 383.00 474.75 7842.25 0.0 0.0 

Enthalpy 14.25 419.00 455.25 7871.50 0.0 0.0 

Estimated Operating Cost for Cooling at $0.10/kWh for Electricity 
(Chiller COP = 3.50) 

Conv (no econ) $10294 for 102939 kWh (Fans - 21693 kWh, chiller = 81246 

kWh) 
Economy $8503 for 85027 kWh (Fans = 21693 kWh, chiller = 63334 

kWh) 
Enthalpy $9295 for 92949 kWh (Fans = 21693 kWh, chiller = 71256 
kWh) 

Estimated Operating Cost for Humidification (20% min) and Heating at $0.10/kWh for 

Electricity 

Humidification $500 for 5002 kWh 
Heating $ for kWh 

Notes: Min space temp occurred on day 345, max space temp occurred on day 147, 
max cooling load occurred on day 165. 



proportional control plan much more attractive. We can conclude that 
TELBECC has great potential for pinpointing significant energy 
reductions and cost savings before a building's HVAC system is 
purchased. 

VI. CONCLUSIONS 

The TELBECC program analyzes more efficiently and quickly than 
any method used heretofore the possible telephone building environ- 
mental energy use and control options. To pinpoint the most econom- 
ical energy-conservation plan, the program analyzes multiple plans at 
minimal cost and with minimal expenditure of time. The program 
calculates the energy consumed every quarter hour by the HVAC in 
regulating the environment under changing weather conditions. It 
computes the required energy from the physical characteristics of the 
building envelope, such as the U factor, internal heat generation, 
geographic location, orientation of the building, and the dry-bulb 
temperature standard. In order to make it feasible to calculate by 
computer, we employ a simplified recursive computation procedure 
using time series. For each of the illustrative problems in Tables II 
through VI, the procedure produced monthly projections; yet it took 
less than 40 seconds to calculate results on an IBM/3033 computer. 
From the examples, we see the advantages and disadvantages of both 
the intermittent and proportional control plans, as well as the signif- 
icant savings obtained from increasing the range of the dual or 
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Table Vlb — Intermittent control of building space air temperature — 
Phoenix, Ariz., case 5 in Table I 

Number of Hours at Specified Relative Humidity (No Humidity Control) 

<10% 10-15% 15-20% 20-55% 55-60% >60% 

Conv(NoEcon) 168.50 611.50 1430.00 6550.00 0\0 0.0 

Economy 199.00 627.00 1432.25 6501.75 0.0 0.0 

Enthalpy 185.50 648.50 1375.00 6551.00 0.0 0.0 

Estimated Operating Cost for Cooling at $0.10/kWh for Electricity 
(Chiller COP = 3.50) 

Conv (no econ) $10527 for 105274 kWh (Fans = 22412 kWh, chiller = 82862 

kWh) 
Economy $8512 for 85123 kWh (Fans = 22412 kWh, chiller = 62711 

kWh) 
Enthalpy $9205 for 92054 kWh (Fans = 22412 kWh, chiller = 69642 

kWh) 

Estimated Operating Cost for Humdification (20% min) and Heating at $0.10/kWh for 

Electricity 

Humidification $1399 for 13991 kWh 
Heating $0 for kWh 

Notes: Min space temp occurred on day 18, max space temp occurred on day 225, 
max cooling load occurred on day 204. 

extended wideband temperature limits for unoccupied times. In the 
larger view, we can understand how TELBECC can significantly 
contribute toward the operating companies' energy-conservation plan 
for future savings. 
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APPENDIX 

Calculating Inside-Wall Temperature Time Series Coefficients bi, a„ a'-, of 

Eq. (1) 

The basis for computing the time-series coefficients is the z trans- 
form, 8 a discrete function transformation. This transformation is 
applied to time functions sampled at regular intervals of time. The z 
transform has the same role in discrete systems that the Laplace 
transform has in continuous systems analysis. 

Let us consider a continuous function of time f(t). When the 
function is sampled at regular intervals A, the output consists of a 
train of pulses, as illustrated in Fig. 1. We defined the z transform of 
this output as a polynomial in powers of z _1 in the following: 

/(0) + /(A)z- 1 + /(2A)z" 2 + /(3A)z- 3 + • • . . (8) 

The successive outputs of the sampler are the coefficients of the 
successive powers of z _1 in the z transform. 

A linear system is characterized when its response to an elementary 
input (such as a pulse, a unit step, or, as will be adopted here, a unit 
ramp) is ascertained. This is nothing more than obtaining a transfer 
function of the system. If both input and output of the system are 
expressed in terms of their z transforms, the ratio of output/input is 
the z transform of the system. If we assume that such a transfer 
function, G(z), can be found and that it can be expressed as the 
quotient of two polynomials in z -1 , then 

rt \ - ^ - °° + ai*" 1 + Q 2 2 " 2 + • • ' a J z ~ J (a) 

U(Z) ~ D(z) " b + fciz" 1 + b 2 z~ 2 + • • • bf ' 

It follows that the z-transform of the output O(z) resulting from an 
arbitrary input I(z) is represented by 

0(z) = G(z)I(z) or (10) 

0{z)D{z) = N(z)I(z). (11) 

Since both sides of (11) are polynomials, the coefficients of the various 
powers of z" 1 must be the same on both sides of the equation. If, say, 
the coefficients of z' n are equated, eq. (11) yields 

b O n = Ooln + aj n -i + a 2 In-2 + ' ' ' + Cljln-j 

- [biO n -i + b 2 O n - 2 + • • • + b p On-p]. (12) 

The subscript n on and J indicates the value of the function at 
t = nA; i.e., 0„ ■ O(nA), the coefficient of z~ n in the z transform of 
0(z). This expression relates the output at any time t = nA to the 
input at that time and the input and output at earlier times. The 
coefficients ao, • • • fly and b , • • • b p contain all the characteristics of 

2956 THE BELL SYSTEM TECHNICAL JOURNAL, DECEMBER 1983 



the system. With the properties of the z transform described above, a 
method for determining the z transform or time-series coefficients for 
the inside building wall temperature follows. 

If we consider the outside building walls and roof structure as 
homogeneous flat slabs (Fig. 2), the temperature in the slab adheres 
to the following equations: 

b 2 T s 8T 8 
* 8x 2 St ' 

T„(x, 0) = 0, 

k ^ (L, t) = -ho[T.(L, t) - To(t)], 
ox 

k ~ (0, t) = h[T s {0, t) - T a (t)l (13) 

where 

T s (x, t) (°F) = temperature in the slab, 
L(ft) = slab thickness, 
fc(Btu/hr — ft — °F) = thermal conductivity, 
k = K/pc(h 2 /hr) = diffusivity, 

pc(Btu/ft 3 ) = volumetric heat capacity, 
/io(Btu/hr — ft 2 — °F) = outside-wall heat transfer coefficient, 
/ii(Btu/hr — ft 2 — °F) = inside-wall heat transfer coefficient, 
To(t) = outside sol-air temperature, and 
T a (t) = inside building air temperature. 

It is convenient to use the Laplace transform 
T s (x, p) = I T s (x, t)e- pt dt 

to eliminate the independent time variable t in eq. (13). Then the 
solution for the inside wall surface (x = 0) temperature in terms of the 
transform parameter p assumes the form: 

^ (Q _ hj[kq cosh(qL) + hp sinh(gL)]T a (p) + hokqT (p) 
hi[kq cosh(qL) + ho sinh(qL)] 

+ kq[kq sinh(qL) + ho cosh(gL)] 

(14) 

where q = (p/*) 1/2 . 

Letting T a (t) and To(t) be unit ramp functions and inverting eq. 
(14) back to the real-time domain by using standard residue theory in 
the complex plane, the solution for T s (0, t), the temperature of the 
inside surface, is expressed as 
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T s (0, t) = T< x) (0, t) + 7?>(0, t), (15) 

where T { s l) (0, t) is the portion of the solution dependent on the outside 
sol-air temperature, and T< 2) (0, t) the part dependent on the building 
space-air temperature. These temperatures are explicitly: 



^"(O, t) = Bo 



_ V_ (3Bj + BpBj + 3B + 6) 

B- + HoB; + So V 6/c (Bi + BoBi + B ) 

(16) 



2L 2 . e -4*» 



k „-i a 2 n \[(Bi + BoB, + B,) - a*] 

• cos a n — a n (2 + Bo + BOsin a n 



} 



and 



, r i a 

' \Bi + BoBi + Bo V 



'/r(0,0 = B i \- , 7> ~ — (^(3 + B ) 



i; (1 + Bq)(3Bj + BjBq + 3B + 6) 

6/c (Bi + BoBi + Bo) ° 

_ 2L^ - [a n cos «„ + Bpsin a n ]e~ At/L * , 

k „=, ^[(Bi + BiB + Bo - a 2 ,] 

•cos a„ — a„[2 + B + BJsin a„ 



J 



where B = -?-, B; = -J-, and a n are roots of the transcendental 

k k 

equation 

a 2 n - BpBj 

COt "* ' «„(B + BiY »-***>'•" 

Equations (15) through (17) contain all the ingredients for forming 
the 2-transform transfer functions for the inside-wall surface temper- 
ature. These are obtained by forming the ratio of output/input z 
transforms as per eq. (9): 

G<»(*)-^^ and G< 2 >(z) = ^J>. (18) 

T (t) and T a (t) were taken as unit ramp functions, and therefore their 
z transform from Ref. 8 is given as 

T (2) = TaU) = 2(1 _ A 2 - 1)2 - (19) 

The sampling interval is A. The use of the input ramp function 
amounts to linear interpolation between the discrete values given by 
the z-transform coefficients. 
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The 2 transforms of both T^iO, z) and T { s 2) (0, z) are similar in form 
and, with the aid of the table of z transforms given in Ref. 8, can be 
expressed as 

4(1,2) R<W>A oo p(W» 

T™(0, z) = f-^ + * * + S -^—t , (20) 

1 — Z 2(1 - 2 T 7=1 1 - SjZ 



where 



^ (1> = , tn . n ° n ■ r, x 2 ( 3S i + fi ofii + 3B + 6), 



-BpL 2 
6/c(Bi + B Bi + Bof 



B d) = "<> 



y3j + BoB; + Bo 

(1) _ _ 2L 2 B / 1 

a?{[(B ; + B B X + Bo) - a]] 
\ -cos otj — a,(2 + B + Bi)sin a,} 

A (2, = B, 



B; + BqBj + Bq 



L * to _, d x L2 d + R o)( 3S i + fi i fi o + 3B + 6) 

— (0 + Do) — — 



,6k 6k (Bi + B Bi + Bo) 

H rJ> 



(2) _ Bj(l + Bq) 



B\ + BoB-, + Bo 

ri2) _ 2BiL 2 I a, cos ctj + Bq sin a, \ 

Cj " " ~T~ ^[(fii + ^Bo+Bo-aJ)] »' an 

V -cos a; — a,(2 + Bo + B;)sin a,- 

Equation (18) can now be expressed in the form of a ratio of polyno- 
mials in z _1 : 

N {h2) (z) 
G {Z) =W) ' (21) 

where, from the results of eqs. (19) and (20), 

N (1 - 2) (z) = (a * Z{1 ~ a Z 1} + B ll A n (1 - sjz- 1 ) 

2(1 - 2" 1 ) 2 °° 
+ U A ) 2 Cff* H (1 - SyS" 1 ), 

D(2) = n (i - sjz- 1 ). 
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Equation (21) is the form of eq. (9); consequently, the 6, coefficients 
that are derived, as in eq. (12), from the coefficients of the polynomial 
D(z) can be generated by a recursive scheme given as 

bo = i, b p = o, b { r 1] = b? - Wi&S n = 1, • • • N (22) 

The number of Sj terms needed to obtain the desired degree of accuracy 
for the bi coefficients is indicated by the index n = N, which in most 
instances should not exceed 20. 

The a, and a[ coefficients in eq. (1) came from N ll) (z) and N {2) (z), 
respectively, by expanding these functions into polynomials in powers 
of 2" 1 ; i.e., 

N ll) (z) = dz" 1 + a 2 z~ 2 + • • • 

N l2) (z) = a[z~ l + a' 2 z- 2 + ■•■ 

The desired coefficients are sorted out. 
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